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Abstract 



This paper deals with the investigation of the solution of an unified fractional reaction-diffusion equation of 
distributed order associated with the Caputo derivatives as the time-derivative and Riesz-Feller fractional 
derivative as the space-derivative. The solution is derived by the application of the joint Laplace and Fourier 
, transforms in compact and closed form in terms of the H-function. The results derived are of general nature 

and include the results investigated earlier by other authors, notably by Mainardi et al. [23,24], for the 
fundamental solution of the space-time fractional equation, including Haubold et al. [13] and Saxena et 
al. [38] for fractional reaction-diffusion equations. The advantage of using the Riesz-Feller derivative lies in 
the fact that the solution of the fractional reaction-diffusion equation, containing this derivative, includes 
the fundamental solution for space-time fractional diffusion, which itself is a generalization of fractional 
diffusion, space-time fraction diffusion, and time-fractional diffusion. These specialized types of diffusion can 
1-^^ ' be interpreted as spatial probability density functions evolving in time and are expressible in terms of the 

, H-function in compact forms. The convergence conditions for the double series occurring in the solutions 

are investigated. It is interesting to observe that the double series comes out to be a special case of the 
Srivastava-Daoust hypergeometric function of two variables given in the Appendix B of this paper. 

Key words: Mittag-Leffier function, Riesz space fractional derivative, Caputo fractional derivative, tele- 
graph equation, Laplace transform, and Fourier transform 
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0> ' 1. Introduction 

o. 

In recent papers, several authors have demonstrated the applications of reaction-diffusion models in pattern 
formation in biology, chemistry, and physics, in this connection, refer to Metzler and Klafter [26], Murray 
[28], Kuramoto [21], Wilhelmsson and Lazzaro [43], Hundsdorfer and Verwer [19], Gafiyuchuk [9,10], Guo 
■ and Xu [12], Chen et al. [2], Engler [5], Pagnini et al. [33], and Huang and Liu [18]. These systems indicate 

, that diffusion can produce the spontaneous formation of spatio-temporal patterns. For details, see the work 

' of Cross and Hohenberg [3], and Nicolis and Prigogine [30]. A general model for reaction-diffusion systems 

is investigated by Henry and Wearne [15,16], Henry et al. [17], Diethelm [4], Saxena et al. [37,38,39,40,41], 
Mathai et al.[25], and Mainardi [22]. 

The object of this paper is to derive the solution of an unified model of reaction-diffusion system (2.1), 
associated with the Caputo derivative as the time-derivative and the Riesz-Feller derivative as the space- 
derivative. This new model provides the extension of the models discussed earlier by authors, including the 
models discussed by Pagnini and Mainardi [33] . Fractional order sub-diffusion is discussed by Naber [29] . 

2. Solution of Unified Fractional Reaction-Diffusion Equations 

Theorem 1. Consider the one- dimensional unified fractional reaction- diffusion equation of distributed order 
aD^N{x, t)+a oD^N{x, t) ^ X ^D^Nix, t) + U{x, t) (2.1) 
xGi?, t>0, 0<a<l, 0<^<1,0<7<2 (2.2) 
with initial conditions 



1 



N{x,0) = f{x), X e R, lim N{x,t) = 0, t > 0, (2.3) 

X— >itoo 

where X is a diffusion constant, o-D*" md qDJ^ are the Caputo derivative operators of orders a and f3 
respectively, xDe^ is the Riesz-Feller derivative of order 7 and skewness 9; \9\ < inin{7, 2 — 7} and U{x,t) 
is a nonlinear function which belongs to the area of reaction- diffusion, then there holds the following formula 
for the solution of (2.1): 

00 I Y poo 

/ \r ft foo 

X (2-4) 

where > 0, 3?(/3) > 0, ^{a - /3) > and 6 = A 'l'^(fc). 

Proof. If we apply the Laplace transform with respect to the time variable t, Fourier transform with 
respect to space variable x and use the initial conditions (2.2), (2.3) and the formula (A3) and (A14), then 
the given equation transforms into the form 

s"7V*(A;, s) - s"-V*(fc) + asl^N^k, s) - as^-'^f*{k) = -A *^(fc)7V*(fc, s) + U*{k, s), 

where, according to the convention followed, the symbol (•) will stand for the Laplace transform with respect 
to the time variable t and * represents the Fourier transform with respect to the space variable x. Solving 
for N*{k, s) yields 

^> sa + asP + X mk) ^ s'^+axff + X ^k) ' ^ ' 

where b = X ipj{k). On taking the inverse Laplace transform of (2.5) and applying the formula (A14), it is 
found that 



N*{k,t)=f*{k) 



5:(-a)'-i(«-«'-i;;+t,),^,(-6t«) 



r=0 -I 

„t 00 

+ / C/*(fc,t-0^(-a)'^r+(«-«'-i£;;+J+(„_^)^(-6r)d6 (2.6) 

The required solution (2.4) is now obtained by taking the inverse Fourier transform of (2.6). This completes 
the proof of Theorem 1 . 

Alternative form of the solution (2.4) 

By using the series representation of the generalized Mittag-Leffler function E^^^{z) in (A15), the expression 

00 

t"-^^(-a)^t("-«'^i;;+V(«-;3).-.+i(-^*")> 

r=0 

can be written as 
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r=0 u=0 



r\u\ r(a — P + 1 + {a — l3)r + au) 



— »-'l:0;0 



-at°'-'^,-bt' 



[1:1;1]:-;- 

[a— p-\-l:a— ^;a]:- 



(2.7) 



where S{-) is the Srivastava-Daoust hypergeometric function of two variables [42]. The definition of this 
function is given in Appendix B. Hence, Theorem 1 can be stated in terms of the Srivastava-Daoust hy- 
pergeometric function of two variables in the following form: Under the conditions of Theorem 1, the 
one- dimensional fractional reaction-diffusion equation 



oD^N{x, t) + a aD^N{x, t) = X ^DlN{x, t) + U{x, t), 
has the solution given by 



(2.8) 



1 f°° 

N{x,t) = —J^ r{k)exp{-ikx) 



cl:0;0 
'-'l:0;0 



-at"-'^,-bt°' 

[1:1;!]:-;- 

— /5;a]: — ; — 

+ ir / / U*{k,t-C)exp{-ikx) 



[1:1;!]:-;- 
[l:a-0;a]:-;- 



dk 



ol:0;0 
^ ''l:0;0 



-af-^^-bt" 



[1:1;!]: 



\a:a—p:a\: 



dk d^, 



(2.9) 



where 3?(a) > 0, ^{/3) > 0, ^{a - ^) > and 6 = A x*^(fc). 



Note 1. By virtue of the Lemma given in Appendix B, the double infinite power series occurring in Theorem 
1 converge for SR(a) > 0, SR(a - ^) > 0. 

3. Special Cases 

When f{x) = S{x), where 6{x) is the Dirac-delta function, we obtain the following 

Corollary 1.1. Consider the one- dimensional fractional reaction- diffusion equation of distributed order 

oD^N{x, t)+a oD^N{x, t) = X a^D^gN{x, t) + U{x, t) (3.1) 
for x £ R, 0<Q!<1, 0</3<l, t > with initial conditions 

N{x,0) =d{x), X € R, lim N{x,t) = 0, t > 0, (3.2) 

where X is a diffusion constant, U{x,t) is a nonlinear function which belongs to the area of reaction- diffusion, 
then there holds the following formula for the fundamental solution of (3.1): 

°° (—nY r°° 

Nix, t)^Y ^ '- / t^''-^> expi-ikx) 

+ ESr^ / C^^"-^^"--' / U*ik,t-OeM-ikx)El+^\^^_^^^{-badk d^, (3.3) 

where 5i(a) > 0, SR(/3) > 0, SR(a - ^) > and 6 = A *^^(fc). 
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If wc set 6 = in (3.1). the Ricsz-Fcllcr derivative reduces to Riesz fractional derivative defined by (A8) 

and it yields the following result given by Saxcna et al. [39]: 



Corollary 1.2. Consider the one- dimensional unified fractional reaction- diffusion equation 



oD^N{x, t) + a oD'^Nix, t) = X ^D^N{x, t) + U{x, t) (3.4) 
for X € R, t > 0, 0<a<l, 0</3<l with initial conditions 

N{x, 0) = f{x), xeR, lim N{x, t)=0, t > 0, (3.5) 

x— >±oo 

where X is a diffusion constant, qD" and o-Of are the Caputo derivative operators of orders a and j3 re- 
spectively, is the Riesz fractional derivative, U {x, t) is nonlinear function which belongs to the area of 
reaction- diffusion, then there holds the following formula for the solution of (3.4)- 



^(^'^) = EV^/ i^"-''^V*(A:)exp(-ite) 

X [E^il-^y+ii-Mkm + a r-'^£;;+i_,)(,^,)^,(-A|fcrr)]dfc 

+E^r^/ r+^^-^-^-w t/*(fc,t-oexp(-zMi?;+i+(„_^)^(-A|fcrr)dfcdc, (3.6) 

r=0 J-00 

where 3?(a) > 0, 3?(/3) > 0, 5R(a - /?) > 0. 

If we further take f{x) = S{x), the above corollary reduces to the following result: 

Corollary 1.3. Consider the one- dimensional unified fractional reaction- diffusion equation 

oD'^N{x, t) + a oD^tN{x, t) = X .xDjN{x, t) + U{x, t) (3.7) 
for X £ R,Q < a < 1, < /3 < 1, t > Q, with initial conditions 

N{x,0) = 6(x), X e R, lim N{x,t) = 0, t > 0, (3.8) 

where X is a diffusion constant, U{x,t) is nonlinear function which belongs to the area of reaction- diffusion, 
then there holds the following formula for the fundamental solution of (3.1). 



N{x, i) = E / exp(-ifca;) 

+ E / U* eM-ikx)E:^-!-a\(a-0)r+ii-Mk\''ndk d^, (3.9) 

where a> (3. 

If we set a = 0, then the theorem reduces to the following result given by Haubold et al. [13, p. 684]. 
Corollary 1.4. Consider the one- dimensional unified fractional reaction- diffusion equation 

oD^N{x, t) = X a:D'^gN{x, t) + U{x, t) (3.10) 
for X € R, t > 0, 0<a<l, 0</3<l, 0<7<2 with initial conditions 

N{x,0) = f{x), x e R, lim N{x,t) = 0, t > 0, (3.11) 
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where X is a diffusion const,ant, and U{x,t) is a nonlinear function which belongs to the area of reaction- 
diffusion, then there holds the following formula for the solution of (3.10). 



1 

N{x, t) = — f*{k)Ec i(-At"*^(A;)) exp{-ikx)dk 

+ 7^/^""'/ U*{k,t-OEaA->^Kik))eM-^kx)dk d^, (3.12) 

where Ea,a{z) is the Mittag-Leffier function. 

Now if we set f{x) = S{x), 7 = 2, ^ = 0, a replaced by 2a and /3 by a, and make use of (A18) then 
the following result is obtained: 

Corollary 1.5. Consider the following one- dimensional reaction- diffusion system 
with the initial conditions 

N{x,0) ^ 6{x), X e R, lim N{x,t) = 0, t > 0, (3.14) 

X— >=too 

where U{x,t) is a nonlinear function belonging to the area of reaction-diffusion. Then for the fundamental 
solution of (3.13), subject to the initial condition (3.14), there holds the formula 

1 f°° 

N{x,t) = , / exp{-ikx)[{a + a)Ea{a t") - (m + a)Ec{fj,t°')]dk 

2-K^{a^ - 46) J-oo 

+ exp{-tkx)U*{k,t-0[Ec.A'^C)-Ea,a{l^C)]dkd^, (3.15) 

where a and fi are the real and distinct roots of the quadratic equation + ay + 6 = 0, given by 

£7= ^(-a+ y(a2^^) andM= ^(-a- y(a2^^), (3.16) 



where b"^ = v'^k'^. 



Next, if we further set U{x,t) = 0, we then obtain the following result which includes many known 
results on fractional telegraph equations, including the one recently given by Orsingher et al. [32] : 

Corollary 1.6. Consider the following one- dimensional reaction- diffusion system 

d'^°'N{x,t) d'^N{x,t) ^d'^N{x,t) 



a- 



'^'^3^'0<a<l (3.17) 



df^ dx'^ 

with the initial conditions 

N{x,0) =6{x), e R, lim N{x,t) = 0, t >0. (3.18) 

>±oo 

Then for the fundamental solution of (3.17), subject to the initial conditions (3.18), there holds the formula 

1 f°° 

N{x,t) = , / exp(-ikx)[(a + a)Ejat°') - (^ + a)EJfJ,t°')]dk, (3.19) 

27rV (a^ - 46) J-oo 

where a and ^ are defined in (3.16), b = u^k"^, and Ea{x) is the Mittag-Leffier function defined in (All). 
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The result (3.19) can be rewritten in the expUcit form as 



where cr and /x are defined in (3.16) and Ea{x) is the Mittag-LefHer function defined in (A17). 

The equation (3.20) represents the sokition of the time-fractional telegraph equation (3.17) subject to 
the initial conditions (3.18), recently solved by Orsingher et al. [32]. It is remarked here that the solution 
given by Orsingher et al. [32] is in terms of the Fourier transform of the solution in the form given below. 
It is observed that the Fourier transform of the solution of the equation (3.17) can be expressed in the form 

4. An Additional Type of Reaction-Diffusion Equation 

Theorem 2. Under the conditions of Theorem 1 with < a < 1 replaced by 1 < a < 2, N{x,0) = f{x) 
and Nt{x) = g{x), and following a similar procedure the solution of the following reaction- diffusion equation 

oD^N{x,t)+aoD^N{x,t) = A ^D^gN{x,t) + U{x,t) (4.1) 

is givn by 

+ J t'^"-<^>g*{k)cxp{-ikx) 



r=0 



fe=0 



X 



d^, (4.2) 

where SR(a) > 0, SR(^) > 0, SR(a - ^) > and 6 = A x*^(fc). 

Corollary 2.1. Under the conditions of Theorem 2 with 9 = the solution of the following reaction- 
diffusion equation 



qD^ + a oD^N{x, t) = X ^DjN{x, t) + U{x, t) (4.3) 



is given by 



ikx) 



+ EV^ / t^"~^^"9*{k)exv{-ikx) 

r=0 



7^0 27r Jo J-oc 
X 0(a-/3).(-^l^r^")dfc dC, (4.4) 
where 3?(a) > 0, SR(^) > 0, SR(a - ^5) > 0. 

Corollary 2.2. Under the conditions of Theorem 2 with < a < 1 replaced by 1 < a < 2 and f{x) = 

g{x) = 5{x), where 5{x) is the Dirac-delta function, the fundamental solution of the following reaction- 
diffusion equation 



is given by 



oD^N{x, t) + a oD^N{x, t) = X a:DjN{x, t) + U{x, t) (4.5) 
N{x,t)=Y^-^ / t("-'^)''exp(-iA;x) 

+ E^^ r U*{k,t-OeM-^kx) 

X K%(o.-0)A-m^k dC, (4.6) 

where Sft(a) > 0, 3ft(^) > 0, ^{a - P) > Q, andb = \ ^^\k). 

As a concluding remark, it is observed that Theorem 2 also holds further if instead of one Riesz-Feller 
derivative, we consider a finite number of Feller derivatives. This result is given in Theorem 3 in the next 
section. 

5. Several Feller Derivatives 

Theorem 3. Consider the following one- dimensional unified reaction- diffusion equation of fractional order 

oD'^N{x,t) + a oD^N{x,t) = [^Xj ^Dj^N{x,t)] + U{x,t), meN (5.1) 

i=i 

for X € R, t > 0, < a < 2, < /3 < 2, < 7j < 2, j = 1, m with initial conditions 

N{x,0) = f{x), Nt{x,0) = g{x), X G R, lim N{x,t) = 0, t > 0, (5.2) 

where Xj > 0, j = l,...,m are diffusion constants, o-D" and o-Of are the Caputo derivative operators of 
orders a and (5 respectively, xDg^,,j = l,...,m are the Riesz-Feller fractional derivatives of orders 7j and 
skewness \9j\ < mini<j<TO[7j, 2 ~ 7j]) j = !> ■■■,'m, respectively. Then there holds the following formula for 
the solution of (5.1): 



N{x,t) = J2^^r^ / t^''-^>f*{k)exp{-ikx) 



X ^;:«+(.-/3).(-^*r)cifc d^, (5.3) 
w/iere 3?(a) > 0, 5R(/3) > 0, 5i(a - /?) > and 6* = 'Z'jli \ "^loi^). 
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Appendices 

Appendix A: Caputo and Riesz-Feller fractional derivatives 

The following fractional derivative of order a > is introduced by Caputo [1[ in the form 



oD^fix, t) = / _i:i_^ dr, m-l<a<m, > 0, m e iV (Al) 



' ' T{m-a)Jo (t-r)«+i-™ 

d"'f{x,t) 



dt" 



ifa = m, {A2) 



where ■^^f[x,t) is the m-th partial derivative of f{x,t) with respect to t. The Laplace transform of the 
Caputo derivative is given by Caputo [1] (also see Podlubny [34], Kilbas et al. [20]) in the form 

m—l 

L{oD^f{x,ty,s} = s^Fix,s)-J2s"-'~'f^'^\x,0+), m - 1< a < m, (AS) 

r=0 

where F{x, s) is the Laplace transform of f{x, t) with respect to t. This derivative is useful in the solutions of 
applied problems connected with anomalous reaction, anomalous diffusion, and anom,a,lous reaction- diffusion 
problems, which are expressible in terms of partial fractional differential equations. In this connection, one 
can refer to the monograph by Podlubny [34], Samko et al. [36], Miller and Ross [27[, Kilbas et al. [20], 
Mainardi [22], Diethelm [4], and recent papers on the subject by Nikolova [31], Naber [29], and Pagnini et 
al. [33]. 

Following Feller [7,8[ it is conventional to define the Riesz-Feller space-fractional derivative of order a 
and skewness in terms of its Fourier transform in the form: 

FUD^f{x);k} = -^i{k)f*{k), (A4) 
where f*{k) denotes the Fourier transform of f{x) with respect to x 

On 

iPl{k) = eyi-p{i{sign k) — ), < a < 2, \e\ < inin{a, 2 - a). (A5) 
Further, when 6 = Q, we have a symmetric operator with respect to x that can be interpreted as 
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This can be formally deduced by writing — (fc)" = — (fc^)"/^. For < a < 2 and \9\ < minjo;, 2 — a}, the 
Riesz-Feller derivative can be shown to possess the following integral representation in the x domain: 

.Dim = ffi±i^l{si„[(„ + .„/21 £' 

+ sm[(a-eW2]j/^"&^l^Mdf). {AT) 

For 9 = 0, the Riesz-Feller fractional derivative becomes the Riesz fractional derivative of order a for 
1 < a <2 defined by analytic continuation in the whole range 0<a<2, aj^l, see Gorenflo and Mainardi 
[11], as 



where 



1 



A 



r—a ^ r2 — a 



2cos(a7r/2)' 

The Weyl fractional integral operators are defined in the monograph by Samko et al. [36] as 



{liN){x) = ^ |_ - C)''-^iV(C)dC, /3 > 



(/^iV)(x) 



{C-xf-^N{QdC, /?>0. 



Note 1. We note that x-Dq ^ pseudo differential operator. In particular we have 



da; 



For 9 = we have 



F{,Dlf{x);k} = ~\k\^f*{k). 
We also need the following result in the analysis that follows: Haubold et al. [13] have shown that 

UI (l,^),(7,|).(l.P)l 



1 



F-'[E^,,{-at^^i{k);x] = —H','i 



(at/3)i/c 



(1 



{AW) 

(All) 
{A12) 

{A13) 



where p= 5R(a) > 0, di{B) > 0, ^{j) > 0. The following results given by Saxena et al. [39] are also 
required in the analysis that follows: From [39[ we have 



„p-i 



s" +as'3 + 6' 



{AU) 



r=0 



where ^{a) > 0, SR(^) > 0, SR(a-/3) > 0, 3?(a-p) > 0, 3?(s) > 0, |^| < 1 and E^^^{z) is the generalized 
Mittag-Leffler function of Prabhakar [35], defined by 



(A15) 



for a, 13, J G C, 3?(^) > 0, SR(7) > 0, where the Pochhammer symbol (a)„ is defined by 

T{a + n) 



{a)n = OL{a + I). ..{a + n-l), (a)o = 1, cc ^ 0, (a)„ = 



T{a) 
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whenver T{a) is defined. When a = 1 in (A15), the expression reduces to the generalized Mittag-Leffler 
function, defined by Wiman [44,45] as 



71=0 



^(13) > 0, »(7) > 0. 



When 7 = 1, (A16) reduces to the Mittag-Leffler function (see, Erdelyi et al.[6]) 



(A16) 



We also have (see, Saxena et al.[39]) 



s2« + as" + b 



- 46) 



[(a + a)£;„(at«) -{fi + a)E„(Mi«)], 



(A17) 



(A18) 



w/iere SR(a) > 0, K (s) > and a and ^ are the real and distinct roots of the quadratic equation +ax+h = 0, 
namely. 



1 



^ = -j^i-O' + \/(a^ -46)) and ji = -(-a - ^ {a^ - 46)). 



(A19) 



Appendix B: Convergence of the double power series 
Lemma. For all a,a,P> 0, there holds the formula 

(l)m+ 



E 

m,n>0 



m!n! (a)am+/3n 



[i:i;i]; 



[a:a;/3];- 



(Bl) 



where S stands for the Sriv astav a- Daoust function (Srivastava et al. [42])- 



The Srivastava-Daoust generalization of the Kampe de Feriet hypergeometric series in two variables is 
defined by the double hypergeometric series as [42, p-151]; 



qA:B-B' 



x,y 



[(a):e,a.]:[(6):*]:[(6'):*' 
[(c):5;6]:[(d):r,]:[(d'):r,'] J 



= E E^™'^ 
771=0 n=0 



mW. 



where 



_ mU r(aj- + m0, + n$,)}{n,ti r(6,- + m^,)}{UU ^b', + n^)} 



(52) 



with the coefficients 9i, ...tOa, ■■■tVD'' > 0. For the sake of brevity (a) is taken to denote the sequence 

of A parameters ai,...,aA with similar interpretations for (6), (d'). Srivastava and Daoust have shown, 
[42, p. 155], that the series (B2) converges for all x,y € C when 



c 

A = l + E'5^- 


D 

^E^i 


A B 

-E^^-E*^>0' 


(B3) 


j=i 
















A' = l + Ee,- 


^Ev 


-E*.-E*/>o- 


m 
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For a, detailed account of the convergence conditions of the double series, the reader is referred to the original 

paper by Srivastava and Daoust [42]. 

Now the left side of (Bl) can be written as 

y^y^ aWn x"^y^ ^r(a)VV r(l + m + n) x^y" 

which establishes (Bl), since for the expression in (Bl), A = q, A' = j3. 
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